Momentum-space calculation of four-boson recombination 
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The system of four identical bosons with large two-boson scattering length is described using 
momentum-space integral equations for the four-particle transition operators. The creation of 
Efimov trimers via ultracold four-boson recombination is studied. The universal behavior of the 
recombination rate is demonstrated. 
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I. INTRODUCTION 

In 1970 V. Efimov predicted that with vanishing two- 
particle binding an infinite number of weakly bound 
states with total orbital angular momentum C = may 
exist in the three-particle system 0, 0] . Such a situation 
takes place if at least two pairs have infinite scattering 
length a — » oo 2]. Away from this limit the number 
of trimers is finite but, depending on a, may be large 
0, Q • This three-particle Efimov effect has an impact on 
the behavior of more complex few-particle systems. In 
particular, for each Efimov trimer the existence of two 
tetramers was predicted in Refs. 0,0. The signatures of 
these states were observed recently in experiments with 
ultracold atoms where the formation of the tetramers led 
to a resonant enhancement of the recombination or re- 
laxation processes 0-0] • 

These first steps in exploring the four-body Efimov 
physics experimentally also call for accurate theoretical 
studies of the four-body systems with large a. Follow- 
ing our recent works devoted to the bosonic atom-trimer 
0] and dimer-dimer 0] scattering, here we focus on the 
four-atom recombination process in the system of four 
identical bosons. The description is based on the exact 
four-particle equations for the transition operators [To| 
derived by Alt, Grassberger, and Sandhas (AGS). For 
the recombination process we need the transition oper- 
ator connecting two- and four-cluster channel states; its 
relation to the two-cluster AGS operators calculated in 
Refs. [8, 9lllllll2|. will be given in the present work. We 
solve the integral AGS equations using momentum-space 
framework. In previous works we demonstrated its reli- 
ability for reactions involving many very weakly bound 
trimers as it happens in the universal regime where the 
accuracy of the coordinate-space methods may be more 
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limited 0, EH, EH- For example, in Ref. 
dieted universal intersections of shallow tetramers with 
the corresponding atom-trimer thresholds while the adi- 
abatic hyperspherical calculations of Refs. 0, 0] were not 
sufficiently precise to find this remarkable feature that 
leads to a resonant behavior in the ultracold atom-trimer 
collisions. Thus, we expect that also in the case of the 
four-boson recombination we will be able to achieve the 
universal limit with higher accuracy compared to the ex- 
isting coordinate-space calculations 0, [l|| ■ We note that 
approximate semi-analytical recombination results have 



been obtained in Ref. [TtJ for a system of three identical 
bosons plus a distinguishable particle. 

An extension of the four-body scattering calculations 
above the four-cluster breakup threshold is a very serious 
theoretical challenge both in atomic and nuclear physics. 
In the momentum-space framework the difficulties arise 
due to complicated singularity structure in the kernel of 
the integral equations. Not all but some of these dif- 
ficulties are present already at the four-cluster breakup 
threshold. In this work we restrict ourselves to a latter 
case which is sufficient to calculate the four-atom recom- 
bination in the ultracold limit. Nevertheless, the present 
work is an important intermediate step towards solving 
the momentum-space four-body scattering equations at 
positive energies. 

We use a system of units where h = 1 and therefore is 
omitted in the equations unless needed for dimensional 
analysis. In Sec. [IT] we recall the AGS equations and de- 
rive the relation between the two- and four-cluster tran- 
sition operators. In Sec. IIIII we present results for the 
four-boson recombination. We summarize in Sec. IIV1 



II. FOUR-BOSON SCATTERING 

We consider a system of four spinless particles with 
Hamiltonian 
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where Hq is the kinetic energy operator for the relative 
motion and the short-range potential acting within 
the pair i; there are six pairs. For the description of the 
scatte ring process in such a system we use the AGS equa- 
tions 10]. They are exact quantum-mechanical equa- 
tions of the Faddeev-Yakubovsky (FY) type. However, 
in contrast to the original FY equations [18| for the wave- 
function components, the AGS equations are formulated 
for the transition operators in the integral form and 
therefore are better suited to be solved in the momentum- 
space framework preferred by us. 

We aim to determine the amplitude for the recombina- 
tion of four free particles into a two-cluster state. Due to 
time reversal symmetry it is equal to the amplitude for 
the four-cluster breakup of the initial two-cluster state. 
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The latter is more directly related to the AGS transition 
operators calculated in our previous works [1, H, EL Ell • 

A. AGS equations 

In a compact notation the AGS equations can be writ- 
ten as 18-component matrix equations [10l.ll4l]. The com- 
ponents are distinguished by the chains of partitions, i.e., 
by the two-cluster partition and by the three-cluster par- 
tition. Obviously, the two-cluster partitions may be of 
3 + 1 or 2 + 2 type; we will denote them by Greek sub- 
scripts. All three-cluster partitions are of 2 + 1 + 1 type 
and are specified by the pair of particles that we will de- 
note by Latin superscripts. For our consideration we need 
the explicit form of the AGS equations for the transition 
operators, i.e., 

U% = (G h Go)" 1 S ap Sji+Y, S*yUi k G t k G U^. (2) 

Here 5 ap = 1 - 5 ap , 

Go = (E + iO-Hq)- 1 (3) 

is the free resolvent with the available four-particle en- 
ergy E, 

U = Vi + v t G U (4) 
is the two-particle transition matrix for the pair i, and 

tff = G - 1 5 jk + h ti G o U; k (5) 

i 

is the subsystem transition operator of 3 + 1 or 2 + 2 type, 
depending on 7. In the summations over the pairs only 
the ones internal to the respective two-cluster partition 
contribute. 

The transition operators U° ap contain full information 
on the scattering process. Their on-shell matrix elements 

Eji<<„<l^pl</>p,J are the amplitudes |T ffp |$ p , n ) 
for the two-cluster reactions. Here 

K,n) ^ GoY^M.n) (6) 

3 

are the Faddeev components of the initial and final chan- 
nel states 

!**»>= I>U>> w 

i 

and n distinguishes between the different channel states 
in the same partition. |$ Pl n) is given by the nth bound 
state wave function in partition p times the plane wave 
with momentum p p ,„ between the clusters, and is nor- 
malized to (®a, n '\$p, n ) = S ap 5 n > n 5(p a ^ - p p , n )- The 
on-shell condition for each channel state relates its bind- 
ing energy, relative two-cluster momentum and reduced 
mass as e Pi „ +p^„/2u p = E. 



The full wave function corresponding to the initial two- 
cluster channel state | Q p , n ) is also determined by the 
AGS transition operators, i.e., 

jjki 

(8) 

The amplitude for the four-cluster breakup reaction 
can be obtained as 

6 

(*o|Top|$,, B ) = (*o|$>i|*„,»>. (9) 

i=l 

The wave function © satisfies also the Schrodinger equa- 
tion, thus, Y?i=i v i\^P,n) = Gq 1 !^,,). Furthermore, 
Gq 1 |$o) = since the four-cluster channel state |$o) is 
an eigenstate of Hq with eigenvalue E. In fact, |$o) is 
a product of three plane waves (each is normalized to 
the Dirac (^-function) corresponding to the relative mo- 
tion of four free particles. Thus, the amplitude for the 
four-cluster breakup of the two-cluster initial state is 

($o|r 0p |$ Pl „) - 5^<* |*i Go c^* Got* G u!;;\f Pjn }. 

jjki 

(10) 

Due to time reversal symmetry it describes also the 
four-particle recombination into a two-cluster state, i.e., 
($ p ,„|T p0 |$o) = ($o|T 0p |$ P ,„). 

If all four particles are identical, there are only two 
distinct two-cluster partitions, one of 3+1 type and one of 
2 + 2 type. We choose those partitions to be ((12)3)4 and 
(12) (34) and denote them in the following by a = 1 and 2, 
respectively. The AGS equations for the symmetrized 
transition operators Up a become 

Uu = P^GotGo) 1 + PziUiGQtGoilii + U^GotGoU^i, 

(11a) 

U21 = (1 + P 3 4)(G iG )- 1 + (1 + P 3i )U x G Q tG Q U xx , 

(lib) 

^12 = (Go<Go) 1 + P3iUiGotGoUi2 + C^GotGoZ-^, 

(11c) 

U22 = (1 + PuPiGotGoUu. (lid) 

Here the two-particle transition matrix t acts within the 
pair (12) and the symmetrized operators for the 3+1 and 
2 + 2 subsystems arc obtained from the integral equations 

U a = P a Go 1 + P a tG U a . (12) 

The employed basis states have to be symmetric under 
exchange of two particles in subsystem (12) for 3 + 1 
partition and in (12) and (34) for 2 + 2 partition. The 
correct symmetry of the four-boson system is ensured 
by the operators P 34 , Pi = P 12 P23 + P13 P23, and P 2 = 
■Pi 3 P24 where P a b is the permutation operator of particles 
a and b. For each two-cluster channel state \& a ,n) = 
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(1 + P a )\<f> a ,n) there is only one independent Faddeev 
amplitude obtained from the integral equation 



\4>a,n) = GotP a \(j> 



(13) 



The symmetrized amplitudes for the two-cluster reac- 
tions are given in Refs. [1, 0, [HI EH- The four-cluster 
breakup amplitude in terms of the symmetrized AGS op- 
erators (fTTj) becomes 

,n) — £*0o ($ |(1 + P 1 ){[l + P 34 (l + Pl)] 
X t Gq JJ\ Gq t Gq I4l a 

+ (1 + P 2 )tGoU 2 G tG U 2a }\<j> 

(14) 

The factors So a depend on the symmetry of the final 
channel state |$o)- For example, Soi = \/6 and S02 = 2 
for nonsymmetrized |$o), while S01 = y/3 and S02 = \/2 
for |$q) that is symmetrized within the pair (12). 

We denote by the four-cluster channel state at 
threshold, where E = and all momenta vanish. Since 
all particles have equal momenta, the channel state |$q) 
is invariant under all permutations, i.e., P a b|$[j) = | <&[]). 
In this particular case the relation (fTfj) simplifies to 



($o|(12t GqU\Gq t G Ui a 

+ 6tG U 2 G tG U 2a )\cl> 



(15) 



We use the momentum-space partial-wave framework 
to solve the AGS equations (jlll) . Two different types of 
basis states \k x k y k z [(l x l y )Jl z ]J r A4) a = \k x k y k z v) a with 
a = 1 and 2 are employed. Here fe x , fc y , and k z denote 
magnitudes of the Jacobi momenta. For a = 1 the Jacobi 
momenta describe the relative motion in the 1+1, 2+1, 
and 3+1 subsystems and are expressed in terms of single 
particle momenta k a as 

1 :(k 2 -k 1 ), (16a) 



[2k 3 -(k 1+ k 2 )], 
[3k 4 - (ki + k 2 + k 3 



(16b) 
(16c) 



while for a = 2 they describe the relative motion in the 
1+1, 1+1, and 2+2 subsystems, i.e., 

1 



(k a -ki), 
(k 4 - ka), 
[(k 4 + k 3 ) 



(ki + k 2 ) 



(17a) 
(17b) 
(17c) 



between all operators in Eqs. (ITT1) . Due to rotational 
symmetry all operators are diagonal in J and indepen- 
dent of M. In addition, U a are diagonal in k z , J and 
! z , P 34 is diagonal in k x and l x , t is diagonal in k y , 
k z and all v 1 and Go is diagonal in all quantum num- 
bers. In this representation the AGS equations for each 
J become a system of coupled integral equations in 
three continuous variables k x , k y , and k z . Such three- 
variable equations have been solved in Refs. [IJ tHI 
for the four-nucleon scattering below the three-cluster 
breakup threshold. The kernel of integral equations fp~T|) 
contains intcgrable singularities arising from each sub- 
system bound state pole of U a ; we isolate these poles 
in different subintervals and treat them by the subtrac- 
tion technique when integrating over k z [11]. However, 
in the energy regime of the four-cluster breakup, E > 0, 
the AGS equations (|TT|) contain additional singularities 
arising from Go, namely, a {k' x k' y k' z v'\Go\k x k y k z v) a — 
5 v>v 5(k' x - k x )S(k' y - k y )S(k' z - k z )/[k 2 x k 2 y kl{E + 1O - 
kl/2/j, ax - k y /2/i ay - kl/2/ia)], where fi ax and fi ay are 
the respective reduced masses. The operator products 
like P 34 P 4 Go render the singularity structure of Eqs. pT|) 
very complicated. It is considerably simpler at E = 
where there is only one singular point k x = k y — k z — 0. 
In this case the integrals involving Go over any of the 
Jacobi momenta kj (except for the trivial integrals in- 
volving momentum (5-functions) are of type 



f(kx, 



iO k^/2fj, ax ky/2fj iQ ,y k^/2fj ia 



(19) 



where the function f(k x ,k y ,k z ) is regular at kj = 0. 
Thus, the the singularity of Go is cancelled by fc| mak- 
ing the integral numerically harmless. Nevertheless, Go 
varies very rapidly near k x = k y = k z = and there- 
fore one should avoid interpolation of Go; this will be 
explained in Sec. Ill Bl 

The four free boson channel state has k x — k y = 
k z = and l x = l y = l z = J = J = 0. Due to the thresh- 
old law the only nonvanishing partial- wave component of 
the breakup amplitude (fTS")) is the one with J = 0, i.e., 



($g|T 0Q |$ a ,„) = 



°|<i> a ,„)/(4^) 2 



(20) 



Although |$q) has only the lj = J = J = compo- 
nent, we emphasize that angular momenta lj and J are 
not conserved. Nonzero lj and J basis states must be 
included when solving the AGS equations (fTTj) . 



The respective orbital angular momenta l x , l v , and l z are 
coupled to the total angular momentum J with the pro- 
jection A4; all discrete quantum numbers are abbreviated 
by v. An explicit form of integral equations is obtained 
by inserting the respective completeness relations 

/>oo 

1 = ^ ^ I \k x k y k z is} a k x dk x k y dk y k z dk za (k x k y k z i'\ 

v JO 

(18) 



B. Separable potential 

The universal properties of the four-boson system must 
be independent of the short-range interaction details. We 
choose the two-boson potential of a separable form acting 
in the S'-wave only, i.e., 



v = \g)H9\8i*o- 



(21) 
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The form factor 

(k x \g) = [l + c 2 (k x /A) 2 ]e^/^ 2 (22) 

is taken over from Ref. [8j. Two very different choices 
of c 2 , namely, c 2 — and c 2 = —9.17, are used in order 
to confirm the universality of the obtained results. The 
strength 

A=-{--[l + |(l+^)l^=r (23) 
irm [a [ 2 V 8 J\ ^/2jt J 

is constrained to reproduce the given value of the scatter- 
ing length a for two particles of mass m. The potential 
(f2"Tj) supports one shallow two-boson bound state at a > 
and no one at a < 0. There are no deeply bound dimers. 
The resulting two-boson transition matrix 

t = \9)r(g\ (24) 

is separable as well with t = (A -1 — (glGolg))^ 1 . This 
allows to reduce the AGS equations (fTTj) to a system of 
two-variable {k y and k z ) integral equations and thereby 
simplifies considerably the numerical calculations. The 
AGS equations with separable interactions are solved in 
the form 

(0|GoHll|01,n) = Pu{9\Pl\<f>l,n) 

+ P34(s|Go[/iGo|5)t( 5 |GoWii|0i,„) 
+ {g\G U 9 G \g)T{g\G(/4n\<t>i,n), 

(25a) 

(g\G U 21 \<f> hn ) = (1 + P 3i ){g\P x \<f>i, n ) 

+ (l + P 3 4){9\GoU l Go\g)r{g\GoUn\cl>i,n), 

(25b) 

(g\G U 12 \h,„) = (ff|l*|&,n) 

+ P34(5|GoC/iGo|5>T( 5 |GoW 12 |0 2:n ) 

+ (<?|GoC/ 2 Go|< ? )r(« ? |GoW22|02,n), 

(25c) 

(.g|G o W 22 |0 2 , n ) = (1 + P 3 4)(5|GoC/iGo|3)T( 5 |G o W 12 |0 2in ). 

(25d) 

We emphasize that with respect to angular momentum 
quantum numbers, the full complexity of the four-body 
problem is still present in Eqs. (|2"5"|) . Basis states up 
to l y ,l Zl J < 3 are included to achieve the partial- wave 
convergence for the solutions of the AGS equations (f2"5f . 
Regarding the continuous variables k y and k z , the treat- 
ment is taken over from Ref. [ill ]. It requires inter- 
polation of (g\GoU a Go\g) in initial and final momenta 
k y . To take care of the Go singularity at kj = 
we factorize (gk y k' z v'\GoU a Go\gk y k z v) = (k' y /2fx ay + 

k' z 2 /2(X a )- 1 U^(k' y , ky, k Z ) {kl^^y + kl/^y 1 5(k' z ~ 

k z )/k 2 z such that V v a u (k' y , k y , k z ) is a very smooth func- 
tion and is interpolated very reliably using spline func- 
tions as in Ref. [llj], while the factors of the type 
(ky/2fj, ay + k z /2/j, a )~ 1 are calculated explicitly wherever 
needed. 



The discretization of integrals in Eqs. (|25|) using Gaus- 
sian quadrature rules leads to a system of linear algebraic 
equations. We reduce its dimension even further by ap- 
plying the so-called /C-matrix technique [3] converting 
the complex equations into the real ones, and by elimi- 
nating the a = 2 component, i.e., by inserting Eq. (I25b[) 
into (I25aj) and Eq. (|25dl) into (|25cl) . To avoid difficulties 
due to possibly slow convergence of iterative methods, 
the resulting linear system is solved using the direct ma- 
trix inversion. Further details of the calculations can be 
found in Refs. [ID, El]. 



III. FOUR-BOSON RECOMBINATION 

The number of the four-atom recombination events in 
a non-degenerated atomic gas per volume and time is 
KiP A where p is the density of atoms and K\ the four- 
atom recombination rate. In the ultracold limit the ki- 
netic energies of initial atoms are much smaller than the 
final two-cluster kinetic or binding energies, and the ini- 
tial momenta are much smaller than the momenta of 
the two resulting clusters, k x ,k y ,k z << p a , n '- Under 
these conditions can be approximated very well by 
the zero-temperature limit K\ calculated at E = with 
k x = k y = k z = in the initial state, i.e., 

Kl = 167r 7 ^/i^ a , n ,|<$°|T ^= |$ a , n ,)| 2 - (26) 

a, to' 

It has contributions from all two-cluster channels (a, n'). 

Concerning the four-atom recombination, the most in- 
teresting regime lies at large negative two-boson scatter- 
ing length a where no shallow dimers exist but the Efimov 
trimers and tetramers cross the zero-energy threshold. 
In our nomenclature we characterize the trimers by one 
integer number n, starting with n = for the ground 
state, and the tetramers by two integers {n,k), where n 
refers to the associated trimer and k — 1 (2) for a deeper 
(shallower) tetramer. The a-dependence of the nth fam- 
ily trimer (b n ) and tetramer (B n _k) binding energies is 
shown in Fig. [1] for a < 0. We denote by a° the specific 
negative value of a where the nth trimer binding energy 
vanishes, i.e., b n = 0. In the universal limit, i.e., for suf- 
ficiently large n, the ratio a° n+1 /a° n = e"/ 80 sa 22.6944 
where s « 1.00624 0. Trimer binding energy 6^ taken 
in the unitary limit l/a = 0, is used to build dimen- 
sionless ratios b n /b™ and Bn^/b™- These ratios become 
independent of the short-range details of the interaction 
provided that its range is small enough compared to the 
size of the few-boson states. This condition is fulfilled 
for high excited states, i.e., for sufficiently large n. We 
checked this independence of the results in Fig. Q] per- 
forming calculations with n = 3 and 4 for two choices of 
the potential form factor (|22p . 

We remind that only two lowest n = tetramers are 
true bound states. All higher tetramers lie above the 
lowest atom-trimer threshold and therefore are unstable 
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FIG. 1. (Color online) Trimer and tetramer binding energies 
as functions of the two-boson scattering length. 

bound states with finite width and lifetime. Their posi- 
tions B nt k shown in Fig. Q] are extracted from the atom- 
trimer scattering calculations at E < as described in 
Refs. [I, [13 . The specific negative values of a — a„ k 
where the respective tetramers emerge at the four free 
particle threshold can be obtained by extrapolating the 
results in Fig. [T] to B n ^ = 0. This procedure yields uni- 
versal ratios 

<i/«! = 0.4254(2), (27a) 
< 2 K = 0.9125(2). (27b) 

The uncertainties, apart from the numerical accuracy, 
arise due to the residual dependence on n > 3 and ci in 
Eq. ([22). 

We will show our results for the four-boson recombi- 
nation rate as functions of the dimensionless ratio a/a^. 
Then in the interval (e~^^ sa , 1) there are exactly n trimer 
states with < n' < n — 1 that contribute to Eq. (121)1) . 
Instead of K\ we build a dimensionless quantity 

Kn = K 4 m/(h\a° n \ 7 ). (28) 

With 1 < n < 4 in Fig. [2] we explore a broad range of the 
two-boson scattering length a and recombination chan- 
nels n' whereas a = 1. The four-boson recombination 
rate varies over many orders of magnitude, but all 
functions of the respective a/a° show qualitatively the 
same behavior. For higher n, i.e., n > 3, there is also a 
very good quantitative agreement, indicating the univer- 
sality of our K n results. Indeed, while n n in Fig. [5] are 
obtained with C2 = in Eq. (|22|) . additional calculations 
with C2 = —9.17 for n — 3 and 4 agree very well with 
the corresponding predictions in Fig. [3J thereby confirm- 
ing our conclusion on n n universality. In contrast, k\ 
shows significant quantitative deviations from the uni- 
versal behavior due to finite-range effects: in the regime 
a/a\ < 1 the only available recombination channel leads 
to the ground state trimer whose size is comparable to the 
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FIG. 2. (Color online) Dimensionless four-boson recombina- 
tion rate K n as a function of two-boson scattering length. 



n 


< ^"n,l/ a n 




a n+l/ a n 





0.4435 


0.8841 


17.752 


1 


0.4412 


0.9162 


21.935 


2 


0.4267 


0.9128 


22.639 


3 


0.4254 


0.9124 


22.691 


4 


0.4254 


0.9125 


22.694 


3 


0.4253 


0.9125 


22.695 


4 


0.4254 


0.9125 


22.694 



TABLE I. Ratios of special two-boson scattering length val- 
ues corresponding to the peaks and cusps of the four-boson 
recombination rate (n > 1) or vanishing tetramer binding en- 
ergy (n = 0). Results are obtained using potential form factor 
with C2 = (top) and C2 = —9.17 (bottom). 

interaction range [l3[ . We note that our non- universal re- 
sults at rt < 1 depend not only on C2 but also on m and 
A; our calculations are performed with m being the mass 
of 4 He atom and A = 0.4 A -1 . 

The four-boson recombination rate has resonant peaks 
at a = a° k where the respective tetramers are at the four 
free atom threshold. The ratios a° fc /a° extracted from 
the peak positions (from B n ^ = in the case n = 0) are 
collected in Table fl] For n > 3 they approach the uni- 
versal limit with high accuracy and are fully consistent 
with the ones in Eqs. (|27[) . Furthermore, the four-boson 
recombination rate has sharp cusps at a = a„, or, equiva- 
lently, at where a new atom-trimer channel 

opens. We list in Table Q] also the ratios a° +1 /a° to prove 
that our numerical predictions converge towards analyt- 
ical result e*/ s ° sa 22.6944 [l|, S] - Furthermore, for n > 3 
we obtain a^y/mb^/h = -1.5077(1). 

At a /a? > 1 not only the ground state but also ex- 
cited trimers are produced via four-atom recombination 
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FIG. 3. (Color online) Relative weights v ni „t for the produc- 
tion of the n'th state trimers via four-boson recombination. 
n = 4 atom-trimer channels contribute in the shown regime. 



Refs. d, HH for the peak positions, i.e., a° x/a% ~ 0.43 
and a° j2 / a n ~ 0.90, are consistent with our results. On 
the other hand, the available experimental data [f| are 
obtained in the system of ultracold 133 Cs atoms with 
deeply bound non-Efimov-like few-atom states and refer 
to n = 0. Therefore it is not surprising that the exper- 
imental results o° 1 /a° l 0.47 and a° j2 / a n ~ 0-84 devi- 
ate from the universal values even more. Nevertheless, 
they remain roughly consistent with the theoretical pre- 
dictions. Thus, our calculations present no improvement 
over the ones of Refs. [J, LL6J when comparing with the 
existing data, but serve as a theoretical benchmark and 
may be valuable for the future experiments performed in 
a truly universal regime. 



IV. SUMMARY 



and contribute to the total rate (pM!) . In Fig. [3] we com- 
pare the relative weights v n ,n' for the production of the 
trimers in the individual states n' when n atom-trimer 
channels arc available in total. The weights are normal- 
ized to Y^n~=o Vn - n ' = 1- The resu lts in Fig. [3] clearly 
demonstrate that more weakly bound trimers are pro- 
duced more efficiently and the four-atom recombination 
process is strongly dominated by the n' = (n — l)th 
channel with the shallowest available trimer. We note 
that a similar conclusion was drawn in Ref. [9] about 
the trimer production in dimer-dimer collisions. Except 
for the vicinity of a/a^ = e _7r / s ° where the shallowest 
trimer disappears, the relative weights v n ,n' depend quite 
weakly on a. Thus, the tetramer states around a = a° k 
enhance the production of all trimers by equal factors. 

As mentioned in the introduction, the recombination of 
four identical bosons has been calculated in Refs. [3, [H| 
using the adiabatic hyperspherical framework. There is a 
good qualitative agreement between our results and those 
of Refs. the shape of the recombination rate, in- 

cluding its resonant peaks and cusps, is reproduced in 
both cases. However, the calculations of Refs. 0, are 
limited to n = 1 where the finite-range effects are not 
entirely negligible and deviations from our universal re- 
sults can be expected much like in the case of our n = 1 
results in Fig. [2] Within this accuracy the predictions of 



We studied ultracold four-boson recombination. Ex- 
act four-particle scattering equations in the AGS ver- 
sion have been used. We related the four-cluster breakup 
and recombination operators to the standard two-cluster 
AGS operators. The scattering equations at the four- 
cluster breakup threshold have been precisely solved in 
the momentum-space framework. The zero-temperature 
limit of the four-boson recombination rate was calculated 
for negative values of the two-boson scattering length. 
We demonstrated its universal behavior in the regime 
with high excited Efimov trimers where the finite-range 
effects are negligible, and determined the positions of res- 
onant recombination peaks caused by the tetramer states. 
We also have shown that most of the Efimov trimers pro- 
duced via the four-boson recombination are in the state 
with weakest binding. Our results are consistent with ex- 
isting experimental data @ and theoretical predictions 
4, 16], but exceed the latter ones with respect to the 
accuracy in the universal limit. In addition, the present 
work constitutes a step towards solving the momentum- 
space four-body scattering equations above the breakup 
threshold. 
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